Some flows in shape optimization 



Pierre Cardaliaguet*and Olivier Ley^ 



Abstract. Geometric flows related to shape optimization problems of Bernoulli 
type are investigated. The evolution law is the sum of a curvature term and a non- 
local term of Hele-Shaw type. We introduce generalized set solutions, the definition 
of which is widely inspired by viscosity solutions. The main result is an inclusion 
preservation principle for generalized solutions. As a consequence, we obtain ex- 
istence, uniqueness and stability of solutions. Asymptotic behavior for the flow is 
discussed: we prove that the solutions converge to a generalized Bernoulli exterior 
free boundary problem. 

Resume. On etudie des hots geometriques lies a des problemes d'optimisation 
de forme du type "probleme de Bernoulli" . La loi d'evolution consideree a la forme 
d'une somme d'un terme de courbure et d'un terme non-local de type Hele-Shaw. 
Notre definition de solution generalisee est fortement inspiree de la notion de solu- 
tions de viscosite. Le resultat central est un principe d'inclusion pour les ensembles 
solutions. Nous en deduisons I'existence, une unicite generique et des proprietes 
de stabilite des solutions. Enfin, nous etudions le comportement asymptotique des 
solutions en montrant qu'elles convergent vers la solution d'un probleme a frontiere 
libre de Bernoulli. 

1 Introduction 

In recent years several works have been devoted to the study of viscosity 
solution for moving boundary problems whose evolution law is governed by 
a nonlocal equation. See in particular [21 [71 El [U [I2l [20l [21]. In this paper, 
we consider subsets of (with N > 2) whose boundary dQ,{t) evolves 
with a normal velocity of the type 
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where A > 0, Ux^^^ is the outward unit normal to (9r2(t) at x, Hx^^'^ is the 
curvature matrix of di}{t) at x (nonpositive for convex sets), F is continuous 
and elhptic, i.e., nondecreasing with respect to the curvature matrix. The 
nonlocal term h is of Hele-Shaw type: 

h{x,Q{t)) = \Du{x)\'^ , 

where u : i^{t) — )■ 1? is the solution to the following p.d.e. 

-An = in n{t)\S, 
u = g on dS, (2) 
n = ond^{t). 

The set S" 7^ is a fixed source with a smooth boundary and g : dS R is 
positive and smooth. We always assume that 5 CC 

The motivation to study such problems comes from several numerical 
works using the "level-set approach" in shape optimization [H |23l |24l [25l 
The idea of these papers is to use formally a gradient method for the 
minimization of an objective function J(Q) where $1 is a subset ofR^. The 
use of the level-set method for building the gradient flow has then the major 
advantage to allow topological changes. Let us underline that this technique 
is up to now purely heuristic. One of the goals of this paper is to justify it 
for some simple shape optimisation problems. 

In order to make our purpose more transparent, a brief description of 
the level-set approach to shape optimization problem is now in order (see 
also the discussion in pLj for a more detailed presentation concerning more 
realistic shape optimization problems) . Consider the problem of minimizing 
the capacity of a set under volume constraints: 

min |cap(f]) with vol(ri) = constant) , (3) 



where 

c.m = / and voKfi) = / 

Jn\s Jn\s 

and u is the solution of ([2]) with instead of il.{t). For any local diffeo- 
morphism 9, we can compute the shape derivatives with respect to 9 of the 
capacity and of the volume. By Hadamard formulas we get 

cap'{n){e)= [ \Du{a)\'^ {9 (a), 1^^) da and vol' (n) (6) = [ {e{a),v^)da . 
Jan Jan 
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Assuming that the optimal shape il. is smooth, the necessary conditions of 
optimahty states that there is a Lagrange multiplier A > such that 



cap'(f])(0) +Avol'(!^)(^) = . 

So it is natural to set 

J^{n) = vol(O) + Acap(n) , 

where A = 1/A. If we choose 6{x) = (—1 + \\Du{x)\'^)i'^ on dQ, then, at 
least formally, we get 

J'^{n){e) = - [ (-1 + X\Du{a)\^ fda < . 
Jan 

Therefore the velocity 9{x) = (—1 + X\Du{x)\'^)i'^ appears as a descent 
direction for the optimization problem ^ and for the set Q. The heuristic 
method for solving ([3]) is now clear: fix an initial position Qq, consider the 
evolution {il.{t))t>o with normal velocity given by ([1]) and F = —1, and 
compute the limit of Q{t) as t — )• +oo: this limit is the natural candidate 
minimizer for ([3|). 

It is worth noticing that problem ([3]) has for necessary condition the 
classical Bernoulli exterior free boundary problem 

Find a set K CC , with S CC K and \Du{x)\ = k for all x G dK, (4) 

where A; > is a fixed constant and u is the solution of ([2]). We refer the 
reader to the survey paper [1^ for a complete description of this problem. 
If one considers a perimeter constraint instead of a volume constraint: 

min jcapfO) with perfO) = constant) , 

one is naturally lead to consider the evolution equation ([T]) with F{h', A) = 
j^^Tr{A) (i.e., the mean curvature). The flow is then formally a descent 
direction for 

Jx{^) = pev{Q) + Acap(r2) . 

Let us underline that this problem has for necessary condition the gener- 
alization of the free boundary problem with curvature dependance (see 

m)- 

Of course all the above computations are only formal: in general, solu- 
tions to the evolution equation do not remain smooth, even when starting 
from smooth initial data. Numerically, this difficulty is overcome by using 
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the level-set approach, which allows to define the solution after the onset of 
singularities. The aim of this paper is to define and study generalized solu- 
tions of the evolution equation, and to investigate the asymptotic behavior 
of the solution as t — )• -|-oo. 

Our concept of solutions is widely inspired by the definition of viscosity 
solution for the mean curvature motion, which corresponds to equation ([1]) 
with F{i>, A) = j^^Tr{A) and A = 0. Motivated by the numerical work 
of Osher and Sethian [22j, a weak notion of solution for this motion was 
introduced in the articles of Chen, Giga and Goto [lO] and Evans and Spruck 
[13| . In this so-called level-set method, the evolution is described as the level 
set of the solution of an auxiliary pde, the level set equation. This equation 
is solved in the sense of viscosity solutions (see [H]). This powerful method 
leads to plenty of results, we refer for instance to the survey book of Giga 
|16j . Note that the level-set approach in shape optimization is a natural-but 
up to now formal-generalization of these ideas. 

As pointed out in [3l [H [26], the generalized solutions obtained by the 
level set approach can also be defined in more geometric and intric ways (see 
also the related notion of barrier solutions introduced by De Giorgi). We 
use here a definition introduced in [2], and used repetitively in [71 13 [9]. In 
the case of the mean curvature motion, Giga [16] proved this definition is 
equivalent to the level-set one. Compared with the already quoted studies 
on viscosity solutions of front propagation problems with nonlocal terms, 
the main novelty of this paper is the fact that we are able to treat signed 
velocities which also involve curvature terms. We learnt recently that a 
similar result (for a Stefan problem) has been obtained by Kim in [21]. 

Our main result is an inclusion principle, which is the equivalent of the 
maximum principle for geometric evolutions. It states that viscosity subso- 
lutions for the flow remain included into viscosity supersolutions, provided 
the initial positions are. For this we have to generalize Ilmanen interposi- 
tion Lemma, which was already the key tool of [TJ |9] . This Lemma allows to 
separate disjoint sets by a smooth (that is C^'^) surface in a clever way. We 
improve this result in two directions (see Theorem 13. 3p . At first we show 
that, when dealing with subsets of R x , the smooth separating hyper- 
surfaces in Rx can be chosen to be smoothly evolving hypersurfaces 
of J?^. Secondly, we build in a carefull way a approximation of these 
evolving hypersurfaces which allows to treat problems with curvature as in 

Let us finally explain how this paper is organized. In Section 2, we de- 
fine the notion of generalized solutions and state the main properties of the 
velocity law. Section 3 is devoted to the interposition Theorems. In Section 
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4 we state and prove the inclusion principle for our generalized solutions. 
As a consequence, we derive results about existence, uniqueness and stabil- 
ity of generalized solutions. Finally, Section 5 is devoted to the asymptotic 
behaviour of the solutions in terms of a generalized Bernoulli exterior free 
boundary problem. 

Acknowledgment. The authors are partially supported by the ACI grant 
JC 1041 "Mouvements d'interface avec termes non-locaux" from the French 
Ministry of Research. 

2 Definitions and preliminary results 

2.1 Definition of the solutions 

Let us first fix sonic notations: throughout the paper | ■ | denotes the cu- 
clidean norm (of R"^ or R^^^, depending on the context) and B{x,R) the 
open ball centered at x and of radius R. If K is a subset of J?^ and x G R^ , 
then dxix) denotes the usual distance from x to K: dxix) = miy^x \y — x\ 
and dx is the signed distance to dK defined by 

dK{x) = I ^^j") , !^ " ^ ^' (5) 
\ -ddK{x) li X e K. ^ ' 

Finally, in the whole paper, if Ki and K2 are subset of R'^ for AT > 1, then 

Ki CC K2 

means that Ki is bounded and, either Ki C mt{K2) or equivalently Ki fl 
W\K^ = $. 

We intend to study the evolution of compact hypersurfaces S(t) = d^{t) 
of R^, where is an open set, evolving with the following law: 

yt>o,xe = hx{x, n{t)) (6) 

where V^^^ is the normal velocity of the evolving set, hx = hx{x, Q.) is given, 
for any set O C R^ with smooth boundary by 

hx(x,n) = Fiu^,H^) + \h(x,n) (7) 

where is the outward unit normal to Q at x, the curvature matrix. 
Throughout this paper we assume that (u, A) e S'^~^ xSn ^ F{v, A) ^ R'ls 
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continuous and elliptic, i.e., nondecreasing with respect to the matrix. Here 
gN-i (denotes the {N — 1)— dimensional unit sphere, and Sn the space of 
A^— dimensional symmetric matrices. Typical examples for F are F{iy, A) = 
— 1 (this corresponds to the flow associated to Bernoulli problem in the 
introduction) or F{u, A) = Tr{A) (for the flow arising in the minimization 
of the capacity under perimeter constraints). As for h, it is a nonlocal 
evolution term of Hele-Shaw type: the example we consider here is 

h{x,n) = \Du{x)f , (8) 

where u : ^ M is the solution of the following p.d.e. 

i) -An = in n\S, 

ii) u = g on dS, (9) 
in) u = on dQ. 

The set ^ / is a fixed source and we always assume above that S CC Q,{t). 
Here and throughout the paper, we suppose that 

i) S C is bounded and equal to the closure of an open set 
< with a boundary, 

ii) g:dS^ (0, +oo) is C^'" (for some a G (0, 1)). 

(10) 

Let us underline that /i(x, Q) is well defined as soon as Q has a "smooth" 
(say for instance C^'°) boundary and that S CC fi. In the sequel, we set 

V = {K^R^ : is bounded and S C int(is:)} , (11) 

where int(ii') denotes the interior of K. 
Prom now on, we consider the graph 

K = {{t,x) eR+ xR^ : xe n{t)} . 

of the evolving sets fl{t). Note that /C is a subset of J?"*" x R^. The set /C 
is our main unknown. We denote by {t, x) an element of such a set, where 
t G R'^ denotes the time and x G R^ denotes the space. We set 

JC{t) = {xeR^ : {t,x) e JC} . 

The closure of the set /C in R^+^ is denoted by /C. The closure of the 
complementary of /C is denoted K: 

K = {R+ X R^) \IC 
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and we set 

= {xeR^ : {t,x) G JC} . 

We use here repetitively the terminology and the notations introduced 
in [HI [9] and [7]: 

• A tube /C is a subset of 1?+ x , such that IC n {[0,t] x R^) is a 
compact subset of R^^^ for any t > 0. 

• A set /C C R^ X 1?^ is left lower semicontinuous if 

Vt > 0, Vx € /C(t), if tra — ^ then 3xn S /C(tn) such that — )• x . 



If s = 1,2 or (1, 1), a regular tube fCr is a tube with a nonempty 
interior and whose boundary has a regularity, and is such that at 
any point (t, x) € dfCr, the outward unit normal '^1^'!^.) = (z^t, i^x) to /C^ 
at (t, x) satisfies 

+ 0. (12) 



is defined by 



• The normal velocity V(j_ of a regular tube /C^ at the point (t, x) G 



where ^^i^-^ = (^^ii i^x) is the outward unit normal to fCr at (t, x). 

A regular tube Kr is externally tangent to a tube JC at (t, x) G /C if 

K. C K-r and (t, x) G dlCr ■ 
It is internally tangent to /C at (t, x) G /C if 

/C-r C /C and (t, x) G . 



• We say that a sequence of C^'^ tubes (/C„) converges to some C^'^ tube 
JC in the C^'^ sense if {JCn) converges to JC and (dJCn) converges to dJC 
for the Hausdorff distance, and if there is an open neighborhood O of 
die such that, if (respectively d^^) is the signed distance ^ to JC 
(respectively to JCn), then (dyc,J and {Ddjc„) converge uniformly to 
dye and Ddjc on O and ||D^d^^||oo are uniformly bounded on O. 
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Remark 2.1 

1. The reason to introduce C^'^ and tubes is clear when looking at ([6]) 
and d?]): h is well defined if ICr is a C^'^ tube (see Section 12. 2p and F is 
well defined if ICr is a tube (to be able to compute the curvature of ICr). 
Therefore, (0) is well defined when ICr is a regular tube and, following 
the ideas of viscosity solutions (see [II]), regular tubes will play the role 
of "test-functions" in the following definition. 

2. For simplicity, we gave all the above definitions with tubes defined for 
all time i > 0. But this is not the point, everything in the sequel is local in 
time, so we use the same definitions with tubes /C where IC{s) is only defined 
in a neighborhood of some fixed t. Note that smooth tubes will always be 
defined locally in time. 

Definition 2.1 Let fC he a tube and Kq £ D be an initial set. 

1. tC is a. viscosity subsolution to the front propagation problem (in short 
FPP) ^ if K, is left lower semicontinuous and fC{t) £ T> for any t, 
and if, for any regular tube ICr externally tangent to IC at some 
point {t,x), with ICr{t) G V and t > 0, we have 

V£) </iA(x,/C,(t)) 

where Vj^^^ is the normal velocity of ICr at (t,x). 

We say that IC is a subsolution to the FPP ^ with initial position Kq 
if IC is a subsolution and z//C(0) C Kq. 

2. IC is a viscosity supersolution to the FPP ^ if K. is left lower semi- 
continuous, and IC(t) C T) for any t, and if, for any regular tube 
ICr internally tangent to IC at some point {t,x), with ICr{t) G T) and 
t > 0, we have 

V(J^) >/jA(x,/C,(t)) . 

We say that IC is a supersolution to the FPP (E^ with initial position 
Kq if IC is a supersolution and if jC{0) C R^\Kq. 

3. Finally, we say that a tube IC is a viscosity solution to the front prop- 
agation problem (with initial position Kq) if IC is a sub- and a super- 
solution to the FPP (with initial position Kq). 

Remark 2.2 The operator h\ defined in d?]) is the sum of a local operator 
F and a nonlocal one h. As in the theory of viscosity solutions, we can local- 
ize arguments related to the local part of the operator. More precisely, 
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regularity of the boundary of the tube is required to compute the curvature 
in F, but only C^'^ regularity is needed to compute the nonlocal part h. 
Therefore, the above definition is equivalent if we replace "for any regu- 
lar tube KLr internally (respectively externally) tangent to K, at some point 
(t, x)..." by "for any C^'^ regular tube /C^ internally (respectively externally) 
tangent to K, at some point (t, x) such that dlCr is in a neighborhood of 
(t,x)..." We will use this equivalent definition in the proof of Theorem 14.11 

2.2 Regularity properties of the velocity h 

We complete this part by recalling the regularity properties of the nonlocal 
term h defined by ([8]) and ([9]) . These results were already given in [7] , so we 
omit the proofs. Here we assume that the set S and the function g satisfy 
assumptions ([TO]) . 

Because of the maximum principle, the function h is nonnegative and 
nondecr easing: if Ki G P and K2 E T) are closed and with a C^'^ boundary, 
if Ki C K2 and if x G dKi n K2, then < h{x, Ki) < h{x, K2). 

Furthermore, h is continuous in the following sense: If Kn and K € T> 
are closed subsets of J?^ with C^'^ boundary such that Kn converge to K 
in the C^'^ sense, if Xn G dKn converge to x S dK, then 



This is a straightforward application of [T71 Theorem 8.33]. 

Next we give a result describing the behaviour of h for large ball: 

Lemma 2.2 For any xq G J?^, there are constants tq > and a > such 
that 



Moreover, the constants rg and a only depend on S and on \\g\\oo- 

The proof is based on standard construction of supersolutions to ([9]) for 
0, = B(0, r), and so we omit it. 

Lemma (j2.2p states that h is small when is a large ball. On the 
contrary, the following lemma means that h is large when "fi is close to 5." 
For all 7 > 0, we introduce 



limh{xn, Kn) = h{x,K). 




S^ = {xeR^ , ds{x) < 7}. 



(14) 



Then, we have 
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Lemma 2.3 There exist 70 > and a constant q > which depends only 
on g and S such that, for all 7 G (0,70), 

h{x,S'y) > — Vx G dS^. 

Proof of Lemma 12.31 Since S has a boundary, we can fix 70 > 
small enough such that ds defined by ([5]) is in S'2'yQ\{d5 < —270}. We 
fix 7 G (0, 7o) and set K = {dg < —7}. We note that dx = + 7 is on 
52'yo\{d5 < —270}. Moreover = 7 on dS and dx = 27 on dS^. Set 

M = max{| A(i|;(x)| | < ds(x) < 7} and m = mm{g{x) \ x G dS} . (15) 

Finally we set Q = S-y and, for /3 = e~^^^^m, we define 

^(r) = /3(e-^^/(4T') - 1) Vr G 1? . 

We claim that 

«(x) = (^(4(x)-(27)2) 

is a subsolution of Q. Indeed, since (^(0) = 0, for all x G dS.y, u{x) = 0. 
From the definition of /3, for all x G dS, u{x) = (y9(— 87^) < m < g. Setting 
rx = d'j^i^) — (27)^, an easy computation gives 

-Au{x) = -^"{rx)\D{dl)\^ - ip'{rx)A{dl). 

But Didj^) = 2dKDdK and {Ddxl = 1. From we get 

-Au{x) < - V(r,)4 + M|/(r,)|. 

A computation of the derivatives of 99 gives 

-An(x) < ^e--^-/(^-^) f 1 - %1 
472 Y 7 

For X G Q\S, we have dxix) > 7 and therefore we obtain —Au{x) < 0. 
Finally n is a subsolution with u > in and n = on dS^/. Thus, for all 
X G dS.y, 

M^;,^^) > |I)n(x)|2 = — . 

QED 

We now recall the main regularity property of the map h: 
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Lemma 2.4 Let R > be some large constant and 7 > 6e sufficiently 
small such that defined by I^14\ ) has a boundary. There is a constant 
9 > 1/7 such that, for any compact set K with C^'^ boundary such that 
Sj C mt{K) and K C i?(0, R — 7), for any v G with \v\ < 1/9 and any 
X G dK, we have 

h{x + v,K + v)>{l-9\v\fh{x,K). (16) 
For the proof, see [TJ Proposition 2.4]. 



3 Interposition theorems 

This part is devoted to interposition theorems in space and in space-time. 
Such results are fondamental in the proof of the inclusion principle. They 
play the same role as Jensen's maximum principle (see [19]) or Ishii's lemma 
(see [m Theorem 8.3]) in the standard theory of viscosity solutions. 

3.1 An interposition theorem in 

Let us start with an interposition result for subsets of . The following 
proposition is a direct consequence of Ilmanen interposition lemma |18) and 
can be found in fP, Proposition 3.7]. 

Proposition 3.1 (Interposition) Let Ki and K2 be two closed subsets of 
R^ , with Ki compact and such that Ki CC K2. Let yi G Ki and 1/2 G 8X2 
be such that 

\yi - y2\ = min \zi - Z2\ . 

zi€Ki,Z2€dK2 

Then there is some open subset Si of R^ with a C^'^ boundary, such that 
El is externally tangent to Ki at yi (i.e., Ki C Si and yi G dT,i) and such 
that S2 := Si + 7/2 — 2/1 is internally tangent to K2 at y2 (i.e., S2 C K2 and 
2/2 G d^2). 

See Figure [T] for an illustration of this proposition. The key point in this 
result is that the smooth set S2 internally tangent to K2 is just a translation 
of the smooth set Si externally tangent to Ki. 

The C^'^ regularity of the sets Si and S2 turns out to be optimal: one 
cannot expect Si and S2 to beC^ in general. Unfortunately the regularity 
will be required in the sequel to be able to deal with curvature terms. In 
order to overcome this difficulty, one can approximate the sets Si and S2 
in the following way: 
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Figure 1: Illustration of the result of Proposition \3.1l 



Theorem 3.2 (Approximation) Let Ki, K2, yi, 1/2, ^1 and S2 be as in 

Proposition l3n\ and 6 > be sufficiently small. Then there exists and 
^2,n open subsets of with C^'^ boundary, converging respectively to Si 
and S2 in the C^'^ sense, there exists yi^n £ Ki and y2,n G dK2 converging 
respectively to yi and y2, and there exists (N — 1) x (A^ — 1) matrices Xi, 
X2 such that 

(i) Tji^n is externally tangent to Ki at yi^n o,nd S2,n is internally tangent 
to K2 at y2,n- 

(a) For i = 1 and 2, is of class C'^ in a neighbourhood of y-i^n with 
limiif3'" — )• Xi, and 

n 

Remark 3.1 

1. Note carefully that the two approximations are not independent because 
of the inequalities (jl7p . which implies in particular that Xi < X2- 

2. By "(5 > sufficiently small", we mean 5 G (0, \yi — ^21/(2 + \yi — 2/2|))- 

The proof of Theorem 13.21 is very similar to the (more difficult) proof of the 
second part of Theorem 13.31 below, so we omit it. 
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3.2 Interposition by regular tubes 

The aim of this part is to extend previous results for subsets T,i,T,2 C 
R X which are regular tubes. The point here is to be able to construct 
tubes satisfying the regularity assumption (fT^ . 

For this we introduce some notations. In Rx we work with the 
norm (where o" > is fixed) 

\{t,x)\„ = (^-^t^+\x\^y . 

For any subset E of R^~^^, we note the distance to E for this norm 
dE{t,x) = inf \{s,y) - . 

For any two subsets Ai and A2 of R^~^^, we define the minimal distance 
between Ai and A2 by 

e{Ai,A2)= inf | (^2, 2:2) - (ti, 2:1)^ • 

{ti,Xi)eAi, {t2,X2)£A2 

We consider the following transversality condition: 

for Ci CC C2 C 1? X 1?^ with Ci compact and C2 closed, 

and for any G Ci and any (82,^2) £ C2, (18) 

if \{si,yi) - {S2,y2)\a = e(Ci,C2), then si > 0,S2 > and yi 7^ 2/2. 

Theorem 3.3 Xei Ci and C2 be such that [T^) holds. Let us fix (si,yi) G 
Ci and (§2,^2) £ C2 wii/i 

l(si,yi) - (S2,y2)|<7 = e(Ci,C2) . 

1. Interposition: There exists aC^'^ regular tube Si, defined on an open 
interval I (see Remark \2. 1[ 2). such that Si is externally tangent to Ci 
at (si,yi), with si £ I, and S2 := Si + (52,^2) — {si,yi) is internally 
tangent to C2 at (52,2/2)- 

2. Joint approximation by tubes: Futhermore, for any (5 > suf- 
ficiently small, there exists C^'^ regular tubes Si^„ and S2,n converging 
respectively to Si and S2 in theC^'^ sense, there exists {si^n,yi,n) G Ci 
and (s2,n) y2,n) £ C2 converging respectively to (si, yi) and (52,^2), and 
there exists {N — 1) x (A'' — 1) matrices Xi, X2 such that 



13 



(i) Tii^n is externally tangent to Ci at {si^n,yi,n) cmd T,2,n is inter- 
nally tangent to C2 at (s2,n;y2,n)- 

(ii) Fori = 1 and 2, is of class C"^ in a neighbourhood of {si^niyi,n) 
with limi??"'"^^'-"^ Xi and 

The proof of this theorem is done in Section 13.51 
Remark 3.2 

1. Inequahty (fT9|) imphes that Xi < X2. Although we only use this latter 
inequality in the sequel, inequality (fTOl) allows to treat equations with F 
depending on x (see for instance [Uj). Let us once again point out that the 
two approximations are not independent because of (I19p . 

2. Thanks to the C^'^ convergence of Si^„ and S2,n to Si and S2 respectively, 
one also has: 

ir = ir = ^S^'"' = ^1?^'^^' 

lim F.^'-" - , = lim y,^'-" - , = _ . = _ _ (21) 

3. By "(5 > sufficiently smah, we mean: 5 £ (0, e(Ci, C^)/(2 + e(Ci, C^))). 



3.3 Existence of the regular interposition tubes 

Let us introduce a new notation: if S is a tube defined on some open interval 
/ (see Remark 12.11 2). then we set 

bd(S) := (J 5S(t) . (22) 

The following result is the key point in the proof of the existence of the 
regular interposition tubes of Theorem 13.31 part (i). 

Proposition 3.4 Let Ci, C2, (si,yi) G Ci, (52,2/2) G C2 he as in Theorem 
\3.3[ There exist a C^'^ regular tube S defined on some interval I and some 
(t,x) G](si,yi), (S2,y2)[ such that 

tel, xe aS(t) and e(Ci , C2) = e(Ci , bd(S)) + e(S , C2) . (23) 
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Above, yi), (s2, y2)[ denotes the open segment joining (si, yi) and (s2, y2)- 



Proof of Proposition 13.41 Let us first fix some notation needed tlirougli- 
out the proof: we set 

• e := e(Ci , C2), 

• E :={(t,x) €R+ xR^ : {t,x) g](si, yi), (s2, y2)[ where 

(si,yi) € Ci and (s2,y2) G C2 satisfy |(si,yi) - (s2,y2)|a = e}, 

• Ap := {{t, x) € X : {t, x) > p and x) > p] for p G (0, 
and, if / is an interval, then Ap{I) = Apf] {I x R^), 

• (t(s),x(s)) := s(si,yi) + (1 - s)(s2,y2) for ah s G (0, 1) and 

(t,x):=(t(l/2),x(l/2)), (24) 

• ■- it - r,t + t) for all r > 0. 

For later use we note that 

d^(i(s),x(s)) = se and <FcAt{s),x{s)) = {l-s)e for ah s G (0, 1) , 

(25) 

because of the definition of (si, yi) and (s2, y2)- Moreover, for a point (t, x) G 
R+xR^, the equality e- (t, x) = (t, x) holds if and only if (t, x) G E. 

We now reduce the construction of the tube S to the construction of a 
suitable function w: 

Lemma 3.5 Let I be a nonempty open interval of R^ and w : Ap(I) — )• R 
be of class C^'^ (for some p G (0, e/2) ) and such that 

e-d''^^{s,y) <w{s,y) <d^cA^,y) V(s, y) G ^p(/) . (26) 

We also assume that there is some 7 G (p, e — p) and some (t, x) G EnAp{I) 
with w{t,x) = 7 and such that 

DMs,y)^0 \/{s,y) e Ap{I) with w{s,y) = j. (27) 

Then the set S = {(s,y) G Ap{I) \ w{s,y) < 7} satisfies the requirements of 
Proposition \3.4\ 
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Proof of Lemma 13.51 Let us first check that S is a tube of class C^'^ in 
the intervall /. Because of assumption ()27p it suffices to show that bd(S) C 
Ap{I) where bd(S) is defined by ([22]). Using ([26]) and the fact that 7 G 
(/9, e — p), we have, for all {s, y) G bd(S), 

d'^^{s,y)>e--f>P and dc^{s,y) > 'y > p . 

Hence {s,y) G Ap{I)- 

Finally we show that ()23p holds. If 'w{s, y) = 7, then (s, y) > 7 while 

(•5) ?/) > e — 7. Hence e(Ci, bd(S)) > 7, e(C2, S) > e — 7 and 

e(Ci,bd(S)) + e(C2,S) >e(Ci,C2) . 

For the reverse inequality let us first recall that e — d'^{t,x) = d^_^{t,x) 
because (t,x) G E. This implies that 7 = (i^^(t,x) > e(Ci,bd(S)) and 

e--f = d^(t,x)> e(C^,S). Hence (El holds. 

C2 

QED 

Next we turn to the construction of a function w satisfying the assump- 
tions of Lemma 13.51 We advice the reader to look at Figure [2] to follow the 
rest of the proof of Proposition 13.41 

The first step is the following result given in [9]: let us set 

Ki := {{t,x) G 1?+ X : d1;^{t,x) < 15e/16} . 
Then, we have 

Lemma 3.6 /9, Lemma 5.1] The function ip{t,x) = dQ^^{t,x) is C^'^ in a 
bounded open neighborhood Oi of E D {Ki\dKi) and 

15e 

e-d'^{t,x)<—-ip{t,x)<d1j(t,x) V(t, x) G 1?+ X J?^ . (28) 

We now show that ip has a nonvanishing spatial gradient in a neigh- 
borhood of E n Oi. Let {t,x) £ E n Oi. Then there exists (si,yi) G Ci, 
(52,2/2) G C2 such that \{si,yi) - {s2,y2)\a = e and {t,x) G](si,yi), (s2,?/2)[- 
Therefore, DLp{t,x) = ((^2,^2) — {siiyi))) From assumption (fT8|) . we 
know that yi 7^ y2- Thus Dxip{t, x) / for all (t, x) G ECiOi. By continuity 
of Dx(p in Oi, there is an open set O2 C Ci which contains E D Ag/Q and 
such that 

r]:= mi-a_\Dxip{t,x)\ > 0. (29) 
(t,2^)eC'2 

We are now going to modify if far away from E. For this we need a 
technical lemma: 
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Lemma 3.7 Let f,g, h be continuous functions in such that g < f < h 
in . Suppose that K := {h = 5} is nan empty and compact and that there 
is some open neighbourhood U of K such that f is C^'^ in U. 

Then, for any fj > and for any open subset U' such that K C U' CC U, 
there is a function ip : R^ — )• R such that: 

(i) V is Cj^l in R^ and is C°° in R\U' , 

(a) g < 4^ < h in R^ and g < i/j < h in R^\U', 

(Hi) \D{ij - f)\<fi inW. 

Proof of Lemma 13.71 Let Ui and U2 be two open subsets of R^ such 
that K C U2 CC Ui CC U' and fix some smooth map 9 : R^ — )• [0, 1] such 
that 61 = 1 in U2 and 61 = in R^\Ui. Then we consider a smooth map 
e,:R^^R such that g < ^ < h in R\U2, 

U - /IIl-(c/'\c/2) < 2\\De\\ "^^ ~ ^/IIl°°{c/'\c/2) < I • 

The construction of such a function ^ is possible because g < f < h outside 
of K and / is C^'^ in U. Then we set 

V'(x) = 9{x)f{x) + (1 - 9{x))^{x) Vx G 1?*^ . 
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Note first that (i) and (ii) obviously hold. As for (iii), it clearly holds in U2 
since ip = f in U2- Moreover, for x G U'\U2, we have 

\D{^P - /)(x)| < \D{^ - f ){x)\ + |/(x) - amD9{x)\ < I + < r? . 

QED 

Next, we apply Lemma O with k := N + 1, U := Oi, U' := O2, f] := r]/2 
where rj is given by (129]) . 



g{t,x) ■.= e-<r^^{t,x)-dX^^{t,x), h{t,x) := dT^it^x) + dX^^{t,x) 

and /(i, x) := 15e/16 - ip{t, x) for {t, x) G J?+ x . We extend /, g and h 
for t < by setting f{t,x) = /(O, x), g{t,x) = g{0,x) and h(t,x) = h{0,x). 
From Lemma 13.61 we have 5 < / < /i in R^~^^. Moreover, from assumption 



([TH]) . the set K := {g = h} = ECiAg/g is compact and contained in (0, +00) x 
and we know from Lemma 13.61 that / is C^'^ in the neighbourhood Oi 
of K. Lemma 13.71 states that there is a map ip : M^'^^ — )• M such that 

' (i) ip is Cl;l in R^+^ and in R^+^\02, 
< iIj < d% in As/a. 

(30) 



(ii) e - < < in A, 



(iii) e-d^ < V<f^&, in A,/8\02, 

(iv) |Z)(V'-/)| <r?/2 inO^. 

Putting together (f29l) and (iv) implies that 



|^xV'(t,a;)| > |I)x/(t,x)| - \D,{^P - f ){t,x)\ > I V(t,x) G O2 . 

We now choose two open subsets Ui and U2 of 1?^ and some r > such 
that 02{t) CC U2 CC Ui (recall that t is defined by (^) and 

|^xV'(i,a:)| > I V(t,x) G X F^. (31) 

We also fix some smooth function 9 : R^ — t- [0, 1] such that 6 = 1 in U2, 
e = Oin R^\Ui and we set 

w{t, x) = e{x)ip{t, x) + (1 - e{x))'tp{t, x) V(t, x) G R^+^ . 

Note that w belongs to Cl;l{R'^+^) nC°°{R x (l?^\[/i)). We claim that we 
can choose r > sufficiently small such that 

r (i) e - < u; < d^c, in^e/4(/r), ^3^) 
1 (ii) iDx'W^I > v/^ in It x Ui. 
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Let us prove the first assertion. Let {t,x) G A^j/^ilr). On the one hand, 
if (t, x) G A^n{I-j-) n {It X [72)5 then 9{x) = 1 and w{t,x) = ip{t,x). Since 
Ag/iilr) C j4e/8, we conclude from pOl fii)). On the other hand, suppose 
{t,x) G Ag/4^{Ir)\{lT X U2). In particular, ^ C'2(i) and (t, x) ^ 

Therefore, from (|3Uf iii)). we have 

Using the uniform continuity of the above distance functions in the compact 
set Ag/^{lT-), we obtain that, for r small enough, 

e-d''^^{t,x) < iP{i,x) < <Fc^{t,x) V(t,x) G A^i^{Ir)\{Ir X U2) . 

Combining with ()30( ii)). we conclude also in this case. 

For the second assertion, we notice that, for G /r x C/i, we have 

\D^{w-i;){t,x)\ < \l-e{x)\\D,{^lJit,x)-iJ{t,x))\ + \i;it,x)-^{t,x)\\D9ix)\ 

with a right-handside smaller than i]/8 provided r is sufficiently small, be- 
cause -0 G Cj^l- Then, for {t,x) G It x Ui, we have, from the choice of Ui 
and T in ([3T]) . 

|i^^i/;(t,x)| > \D^i;{t,x)\ - \D^{w - i;){t,x)\ >r//8, 

which proves the second statement. 

We now fix (7 G (0, e/4) such that 

{t{s),x{s)) £lrX 02 (t) Vs G (1/2- (7,1/2 + . 

This is possible because (t,x) = (t(l/2), x(l/2)) belongs to 02{t). Prom 
(j25|) . an easy calculation gives 7i'(t(s), x(s)) = (1 — s)e. Since w is smooth 
in 1? X (l?^\?7i), Sard Lemma states that we can find a level 7 G ((1/2 — 
a)e, (1/2 + (T)e) such that 7 is a non critical value of in 12 x (J?^\C/i). 
We claim that w and 7 satisfy the requirements of Lemma 13.51 Note first 
that, for s = (e — 7)/e, the point (t(s),x(s)) belongs to £' n Ag/4^{Ir) and 
satisfies w(t(s),x(s)) = 7. Moreover, (p6|l holds from (|32r i)). Finally we 
show that (f27|) holds. Indeed, if w{t,x) = 7 for some (t, x) G A^j^^It)-, then 
either x G t/i, in which case Dxw{t,x) / thanks to ([32lf ii)). or x ^ C/i 
and Dw{t,x) = {0, Dxw{,tx)) 7^ because 7 is a non critical value of w in 
Ir X (12^\c7i). In each case we have Dxw{t,x) ^ 0. This completes the 
proof of Proposition 13.41 

QED 
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3.4 regularization of a C^'^ tangent surface near contact 
points 

The aim of this section is to show the fohowing fact: if a C^'^ surface S 
is externally tangent to a set K at a point y, then it is possible to find a 
C^'^ surface S which is close to S (in the C^'^ sense) and is still externally 
tangent to IT at a point y close to y. Moreover, S is more regular than S, 
namely is in a neighborhood of y. In particular, we can use S as a test 
set to estimate the curvature (see Remark I2.2p . 
Let us give the exact assumptions: 

Let K he a subset of -R'^ for k > 1 and S be an open set with a 
C^'^ boundary which is externally tangent to K at some point 
y £ dK. Let x ^ K he such that y is the unique projection of x onto 
K and p := Ddxix) is the outward normal to S at y. Suppose that (33) 
there is a sequence of points where is twice differentiable 

with first and second derivative denoted respectively pn and Xn , and 
finally assume that pn converges to p while Xn converges to some X. 

Note that, by usual properties of the distance function at differentiability 
points, the projection of x„ onto K is unique and converges to y. We denote 
by yn this projection. 

Proposition 3.8 Under Assumption ^33\) we can find a sequence of open 
sets S„ with C^'^ boundary such that 

(i) Tin is externally tangent to K at yn, 

(a) S„ has a boundary in a neighbourhood of yn, with normal Pn and 
curvature equal to the restriction to (pn)"*" of —{Xn — ^^k), 

(Hi) S„ converges to S in the C^'^ sense. 

Before starting the proof of the proposition, we need two lemmas. The first 
one builds, from the derivatives of the distance function at a point a, a map 
(f) which has a local maximum on K at the point b, projection of a onto K: 

Lemma 3.9 Suppose that a ^ K and that dx is twice differentiable at a. 
Let b be the projection of a onto K. Then, for any a > 0, the (smooth) 
function 

(t>{z) = {DdKia), z-b) + ^{{D^dKia) - alk){z -b),z- b) 
has a strict local maximum at b on K. 
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Proof of Lemma 13.91 For any z & K, we have 

(f){z) = dxiz + a — b) — —\z — b\'^ — dxia) — \z — h\^^{z — h) 

because is has a second order Taylor expansion at a. But, since z ^ K, 
we have dxiz + a — h) < \{z -\- a — b) — z\ = \a — b\ = dxia). Therefore 

Hz) < -^\z - b\^ - \z - b\h{z - b) 

which is negative as soon as z G is sufficiently close to b and z ^ b. 

QED 

From now on, we fix a smooth function 9 : 12+ — t- [0, 1] such that 9 is 
nonincreasing, 9 = 1 on [0, 1/2] and 9 = on [1, +oo). 

We will use several times below the following interpolation. The proof 
relies on straightforward computations, so we skip it. 

Lemma 3.10 Let cj) and Tp some C^'^ functions in some open set O. Let 
y £ O be such that Hy) = ip{y), and let us set, for any p > 0, 

^z) = <P{z)9p{z) + ^{z){l - 9p{z)) where 9p{z) = 9 

Then, for any p> such that B{y,p) CC O, we have 
Up - V'lloo < C{t]p + (Ml + M2)p2) , \\D^^ - Z^Vlloo < C{7] + (Ml + M2)p) 
and 

P^lloo < -{v + {Mi + M2)p) 
P 

for some constant C = C{k) > 0, where we have set r] = \D(l){y) — D'ijj{y)\, 
Ml = \\Dm\oo and Ms = P^V^IU- 

Remark 3.3 

1. The norms || • ||oo are of course taken on ^(y,p), since ^p — 'ip outside. 

2. The key point is that coincides with in a small neighbourhood of y, 
but is not too far from in the full set O provided p and r] are small. 

We are now ready to prove the proposition. 
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Proof of Proposition 13.81 From Lemma 13.91 the function (pn defined by 

(pniz) = {DdK{Xn), Z - yn) + \{{D'^dK{Xn) " -Ik){z - Vn), Z - Vn) 

1 n 

has a strict local maximum at y„ on K. Since y„ is the unique projection of 
Xn onto K, the map V'n(-z) = dxixn) — \z — x„| has a global strict maximum 
on K at Un. Hence 



Cn(^) = (l)n{z)On{z) + llJn{z){l - 6n{z)) where Oniz) = 6 



Vn? 



has a global strict maximum at yn on K, provided we choose p„ > 0, p„, — )■ 
and n large enough. Since (pn is globally smooth with uniformly bounded 
second order derivative, and since ^„ is smooth with uniformly bounded 
second order derivative ouside of the ball B{x,dKix)/2), and since finally 
tpniyn) = (pniVn) = and D^niVn) = DipniVn) = Vn, Lemma ETJO] states 
that C,n and DQn uniformly converge to the function z — )■ dxix) — \z — x\ and 
its derivative respectively, in the set R^\B{x,dK{x)/2), and that ||-D^Cn||oo 
is uniformly bounded in R^\B{x,dK{x)/2). 

Let us now denote by ds the signed distance to S (see ([5]) for a defini- 
tion). Since 9S is a C^'^ manifold, we can find some open neighbourhood O 
of as such that ds is C^'^ in O, with HD^dsHoo bounded in O. For z G R'', 
we define 

d'^iz) = dj:{z) - I3n\y - where f3n > 0, f3n ^ 0. 

Notice that is C^'^ in O, that c?" and its derivative converge locally uni- 
formly to ds whereas jjl^^d^lloo is bounded in O. The advantage of in- 
troducing is that {d^ < 0} is still externally tangent to K at y with 
dK D d{d"' < 0} = {y} (instead, S can touch K at many points). We claim 
that, if we choose /3„ = 2|y — then, at least for n large, 

d"(z) < d"(y„) for any z G K\B{yn, /3„/2). (34) 

Indeed, for z G K\B{yn, 

dr{z)-dr{yn) < -Pn\y- Z\'^ -d^{yn) + lin\y-yn? 

< -/3n(/3n/2 -\y - yn\f + \y- yn\ + Pn\y - yn? 

< -|y-yn|(i-4|y-y„|§), 

which is nonpositive for large n since yn y- 
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We introduce the maps 



Uz) = C„(2)0n(^)+K(^)-d"(yn)](l-en(^)) where dniz) 
and 



an = max{|p„ - /5n/\/2} . 

Let us notice that o"n — t- because y andp„ and D(F'{yn) = DdY,{yn) — 
2fin{Vn — y) converge both to p = Ddj]{y) since Dd^ is continuous at y. 
We now use Lemma [3.101 with QniVn) = [d^{yn) — d^iyn)] = and rjn ■= 
\DCniyn) - D(P{yn)\ = \pn - Ddj:{yn) + 2/3„(y„ - y)| 0. It states that 
and D^^n uniformly converge to and Dds, respectively. Moreover, since 
Vn < CTni the second order derivative of ^„ is uniformly bounded. 

Let us finally prove that the set S„ = < 0} satisfies our requirements. 
What we already proved on ^„ shows that S„ converges to S (in the C^'^ 
sense). From its construction, S„ is smooth in a neighborhood of y„, with 
normal at the point yn equal to Pn and curvature equal to the restriction to 

(Pn)^ of - i/fc). 

It remains to check that S„ is externally tangent to K at yn- It suffices 
to prove that ^n(-z) < for any z £ K, because S,n{yn) = 0. Let z £ K. If 
\z — yn\ < CFn/V^, then ^n{z) = Cn{z) < from the construction of Cn- If 
\z — yn\ > an/V^, then \z — yn\ > Pn/2 and thus, from < 
Since moreover ^„ has a global maximum on K at y„, we finally have 

Uz) = Qn{zyen{z) + [d^{z)-r{ynm-en{z)) 

< UynVeniz) + [d^{yn) " d^iyn^l " On{z)) = . 

In conclusion we have proved that has a global maximum on K at y„, 
and the proof is complete. 

QED 



3.5 Proof of Theorem ISTSl 

The first part of the theorem is an immediate consequence of Proposition 
13.41 : we set Si = S — (t, x) + (si, yi) and S2 = S — (t, x) + (s2, 1/2) where S 
and (t,x) are given by Proposition 13.41 and we check that Si and S2 enjoy 
the desired properties. 

Without loss of generality, we assume that 5 G (0,e/(2 + e)). Let us 
introduce, for all (n, zi, r2, Z2) G 1?+ x 1?^ x 1?+ x K^, 

/(ti,2:i,T2,Z2) = dc^{Ti,zi) + d'^{T2,Z2) + ;^|(ri,zi) - (t2, 2:2)1^ • 
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Then d^^, d"^ and / are semi-concave functions in (1?+ x l?^)^\(Ci U C2)^. 
We claim that / has a minimum at (ti, ^i, f2, ^2), where 

in,zi) = ^(1 + 6){si,yi) + ^(1 - <5)(s2,y2) 

and 

{f2,Z2) = -{l-S){suyi) + -{l + 6){s2,m). 

Indeed, on the one hand, an easy computation shows that f{fi, ^1, T2, Z2) = 
e - 6/2. On the other hand, for ah (n, zi, T2, 22) e x R'^ x R+ x R'^ , 
we have 

f{Tl,Zi,T2,Z2) 

= dci(n,^i) + (i^(r2,^2) + ^|(ri,zi) -(t2,Z2)|' 

= \{n,Zi) - {ti,Xi)\^ + |(r2,22) - {t2,X2)\a + ^K"^!'^!) ~ ("^2,^:2)1^ 
> \{ti,Xi) - {t2,X2)\a - \ {ri,Zl) - {T2,Z2)\a + ^K"^!'^!) ~ ('^2,^:2)1^) 

where {ti,xi) G Ci and {t2,X2) G C2 are the points for which the distances 
d?, (ti,zi) and d'^(T2,Z2) are achieved. It fohows that /(ti, zi, r2, 2:2) > 

e — 6/2 since | (ti, a;i) — (t2, a^2)|(T > e and since, for all r > 0, — r + r^/ (26) > 
—6/2. Finally, (ti, zi, f2, ^2) is a minimum for /. 

Since the semi-concave function / has a minimum at (ti, zi, f2, Z2), Jensen 
maximum principle [19] (see also [I4j) states that one can find a sequence 
of points (Ti,n; ^i,n; T2,n, 22,n) which Converges to {fi, Zi,f2, Z2) and a non- 
negative symmetric matrix A G 52Af+2 such that the functions , dg^ and 

/ are twice differentiable at (ri,„,zi,„), (f2,„,22,„) and (fi,n, zi,n, T2,n, ^2,n) 
respectively and such that 

-D/(ti,„, r2,n, ^2,n) and -D^/(fi,„, f2,n, Z2,„) ^ > 0. 

(35) 

In particular, since Dd'^_^{Ti^n, zi^n) and Dd'^{T2-n, Z2-n) exist, the projec- 
tions of (ti^rn^i.n) onto Ci and C2 respectively are unique, and equal to 
some {si^n-,yi,n) and {s2,n,y2,n)- Note that (si,yi) is the unique projection 
onto Ci of (t, x), and therefore (si,n) yi,n) converges to (si, yi). For the same 
reason, {s2,n-,V2,n) converges to {s2,y2)- Since d% and dS;- are semi-concave 
functions, the matrices D'^d% (ti^„,zi^„) and -D^(i2^(f2,n, ^2,n) are bounded 
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from above: namely (see for instance [6^ Proposition 22.2]) 

^L^Ci(n,n,^i,n) < ,^ _ J n (36) 



and 



DljUf2,n, Z2,n) < - J n ■ (37) 



Using (j35|) . we get 

C2 



It n2^a/- - N^ir ^^>0, (38) 



where the matrix A G S2N is the restriction to 1? x R oi A. In particular, 

if we set ^ = f I , then 

V ^3 A2 J' 

2 

^Lc^Ci (n,n, ^l,n) + I^L^^gr (r2,n, ^2,„) + ^/^V ^ ^1 + ^2 > 

and therefore D^^d^^ (ri_„, zi^„) and D'^^d'^{f2^n, Z2,n) are in fact bounded. 
So, after relabelling all the sequences, we can assume that the restriction 
of —D'^^dQ^{fi^n, zi^n) to {Dxd'Q^{fi^n, zi^n))'^ Converges to some matrix Xi 
while the restriction of D'^^d'^^ {f2 n, Z2 n) to (Dxd'^^(f2 n, Z2 n))'^ converges 



to some X2. Note that, from (]38p we have 
Moreover, since 



dcS'^l,n, Zl,n) dc^{Ti,Zi) = ^ 

and 

C?gr(r2,n,^2,n) ^ C^g^(T2,^2) = ; 



we get from ([361 [3j 

V X2 j - 1(137)^2(^-1) < f2iN-l) 

because (5 < e/(e + 2). So we have proved (fT9ll . 
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We now apply Proposition 13.81 to the sets Ci and Si. Assumption (p3]) 
holds since the set Si is externally tangent to Ci at (si,yi). Moreover, the 
point (si,yi) is the unique projection of the point (ri,zi) onto Ci. The 
points (Ti,n,^i,n) are points of twice differentiability of d^^, converge to 
(fi,zi) and have a unique projection (si^nji/i.n) onto Ci. Therefore we can 
find a sequence of sets Si^„ with C^'^ boundary, such that Si^„ is externally 
tangent to Ci at {si^n,yi,n) and has a boundary near (si^„,yi^„). Note 
that, since Si is a C^'^ regular tube and since the sets Si^„ converge to Si 
in the C^'^ sense. Si „ are also C^'^ regular tubes provided n is sufficiently 
large. 

From Proposition 13.8( 11). we also have that the curvature matrix of Si „ 
at {si^n,yi,n) is equal to the restriction of — (D^d^^ (ri^„, zi^„) + ^I^^i) to 
the tangent space of Si^„ at (si,n, yi,n)- Let us denote by Xi^n the restriction 
of this curvature matrix to . We notice that Xi.„ converges to Xi. 

In the same way, applying Proposition 13.81 to the complementary of the 
tube S2 which is externally tangent to C2 at (52,2/2), we can find a se- 
quence of C^'^ tubes S2,n which are externally tangent to C2 at some points 
{s2,n,y2,n), such that S2,n are of class near {s2,n,y2,n), and such that the 
curvature matrix X2^n to S2,n(^2,n) at y2,n converges to X2- 

QED 

4 The inclusion principle 

4.1 Statement of the main theorem. Existence, uniqueness 
and stabihty 

Theorem 4.1 (Inclusion principle) Let < Xi < X2 be fixed, fCi be a 
subsolution of the FPP ^ with speed hx^ on the time interval [0,T) for 
some T > and IC2 be a supersolution on [0,T) with speed hx^. If 

^(o)n£^(o) = 0, 

then, for all t G [0,r), 

^(t)n^(t) = 0. 

Before proving Theorem 14. H we recall some applications of such inclusion 
principle, omiting the proofs which are easy adaptations of those of and 

Concerning the existence of solutions, we have the following 
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Proposition 4.2 For any initial position Kq, with S C int(i^o) cLnd Kq 
hounded, there is (at least) one solution to the FPP ^ for h\. 

Moreover, there is a largest solution and a smallest solution to this prob- 
lem. The largest solution has a closed graph while the smallest solution has 
an open graph in x . The largest solution contains all the suhsolu- 
tions of the FPP ^ with initial condition Kq, while the smallest solution 
is contained in any supers olution. 

In general, one cannot expect to have a unique solution, i.e., the closure 
of the minimal solution is not necessarily equal to the maximal solution. 
However, uniqueness is generic: 

Proposition 4.3 Let {Kx)\yQ be a strictly increasing family of initial po- 
sitions (i.e., Kxi CC K\ for < A' < such that K\ € T) for all A > 0. 
Then the solution of the FPP ^ for h\ with initial position K\ is unique 
hut for a countable subset of the A 's. 

Stability of solutions is expressed by means of Kuratowski upperlimit of 
sets. Let us recall that, if {An)n&i is a sequence of sets in 1?*^, then the 
Kuratowski upperlimit A* = Lim sup„ is the subset of all accumulation 
points of somes sequences of points in namely 

A* := {z G 1?^^ : 3{nk)k(zfq increasing sequence of integers, 3(zfc)fcgN) 



We define A^ as the complementary of Limsup„ j4„. 

Proposition 4.4 // /C„ is a sequence of subsolutions for hx, locally uni- 
formly bounded with respect to t, then the Kuratowski upperlimit IC* = 
Limsup„/C„ is also a subsolution for hx. 

In a similar way, if fCn is a sequence of supersolutions for hx, locally 
uniformly bounded with respect to t, then /C* is also a supersolution for hx- 

4.2 Proof of Theorem 14.11 

Without loss of generality, we assume that /Ci has a closed graph. We argue 
by contradiction, assuming there exists < T* < T such that 



Zk G ^rifc and z = limyfc}. 



k 



(39) 




(40) 
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For a > and e > 0, we consider 

Pa{t) = min d'^{t,x) 

and we set 

T'^" = mi{t > : Pa{t)< ee"* } . 

Recall that the notations I • L and df^ were introduced at the beginning 

of subsection 13.21 Let r > sufficiently small such that S + rB has a 
boundary and S + rB CKZ )Ci{t) and S + rB CC K.2{t) for ah t G [0, T]. We 
also fix i? > sufficiently large such that 

sup \x\ + sup |x| < R — r . 

{t,x)e!Ci, t<T {t,x)£K.2, t<T 

We denote by 6 the constant defined in Proposition 12.41 for R and r. Let 
us recall that 9 > 1/r and that, for any compact set K with C^'^ boundary 
such that Sr C int(i^) and K C B{0, R-r), for any v e with \v\ < 1/6* 
and any x G dK, we have 

h{x + v,K + v)>{l-9\v\fh{x,K). (41) 

We refer the reader to Figure [3] for an illustration of the proof. 











, 












''"'^'i'iA-'^C' 
















T'-" = si 



Figure 3: Illustration of the proof of Theorem \4.1[ 

Lemma 4.5 We can choose e > and cr > sufficiently small so that 

(i) A2(l - > Ai and ea < 1 , 

(ii) T'^" > e, 
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(ill) p^iT^n = ee-^"\ _ 

(iv) for any yi G ICiiT^'"^) and any (52,2/2) G ^2 such that 

pAT'n= min d^(r^'^y) = |(r^-^yi)-(s2,y2)U, (42) 
we /laue yi / 7/2 o'^c^ 52 > 0. 

The proof of the lemma is postponed to Section 14.31 

Prom now on we fix e and a as in Lemma 14.51 and we set si = T'^'^. Let 
us also set 

JC];" = {{s,y) e R+ X : 

3(r, z) G /Ci with r < si and \{s,y) - (r, z)|^ < e(e~^ - e"""^)}. 

Recall that, for any two subsets Ai and A2 of we define the minimal 

distance between Ai and ^42 by 

e{Ai,A2)= inf | (^2, 2^2) - (ti, 2;i)|a • 

(il,xi)GAi, (t2,a;2)GA2 



Lemma 4.6 The set IC^" is closed, with ICl'^{si) = /Ci(si) and 

e{}C'{%}C2) = ee-'\ 

Moreover there exist yi G /C^''^(si) and (^2,2/2) £ K^2, such that 

l(si,yi) - {S2,y2)\a = eilC^" (43) 
and, for such points yi and (^2,^2); we have yi ^ 1/2- 

The proof is postponed. Lemma 14.61 is a kind of refinement of Lemma 
14.51 (iv). Note that the proof of Lemma [4. 51 and Lemma [4. 61 only use the fact 
that /Ci(0) n Al2(0) = and that /Ci and IC2 are left lower semicontinuous. 

Next we give an estimate of the normal velocity of the tube tCi in terms 
of the normal velocity of fC^": 

Lemma 4.7 Assume that a C"^'^ tube E is externally tangent to Kf^" at some 
point {s,y) G dlC^'^ . Let {f,z) G JCi be such that 

f<si and \{s,y) - {f,z)\„ < e{e-^ - e-''). 
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Then there is a C ' tube S which is externally tangent to the tube /Ci n 
([0, f] X K^) at (t, z) and such that 

^(|,-)>T/(J5)+ee-«^ and S(f ) = S(s-) + (z - y). 

If furthermore S is of class in a neighbourhood of {s,y), then T, is also 
in a neighbourhood of{f,z). 

We postpone the proof. We are now ready to use the interposition Theorem 
13.31 that we apply to Ci := /C^'"^ and C2 ■= K,2- Note that condition (fTSj) 
holds thanks to Lemma WM Let us fix G IC^" and (52,^2) £ ^2 with 

|(si,i/i) - (S2,y2)| = e{lC{",lC2) ■ 

From Theorem 13.31 we know that there exists a regular tube Si with C^'^ 
boundary such that Si is externally tangent to /C^''^ at (si,yi) and S2 := 
5^1 + (•52)2/2) — (•51, yi) is internally tangent to IC2 at (•§2,^2) (see Figure [3]). 

Futhermore, there exist C^'^ regular tubes Si^„ and S2,n converging re- 
spectively to Si and S2 in the C ' sense, there exist (si,n,yi,n ) G 1C{^ and 
(s2,n,y2,n) £ ^2 Converging respectively to (si,yi) and (s2,?/2), and there 
exist (A^ — 1) X [N — 1) matrices Xi, X2, such that 

(i) Si^n is externally tangent to IC^" at (si,n,2/i,n) and S2,n is internally 
tangent to IC2 at (s2,n,y2,n), 

(ii) For i = 1,2, Sj^„ is of class in a neighbourhood of {si^n-,yi,n) with 

lim^.£:(^^-") = hmz^gr^--") = vl^^^^) = v^-(^-\ (44) 

limy.^''" ,=limy,^''" , = VT^ ^=yT^ ^, (45) 
and, for i = 1, 2, 

lim fI;';':^, „ ^ Xi, with Xi - X2 < . (46) 

n 

Since /C2 is a supersolution for h\^ and S2,n is internally tangent to /C2 
at (•S2,n,y2,n), we have 

KZm..) ^ + A2My2,n,S2,n(s2,n)). 
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Letting n — t- +oo gives, from the continuity property of h, from 
and (06]) 

^Sly2) ^ n^m^''^^X2) + X2h{y2,Ms2)) (47) 

We now establish a symmetric inequality for Si. For any n, let (ti „, „) G 
fCi be such that 

Ti,n < Si and - {n^n, zi,n)\a < e(e"^i'" - e'"^). 

Note that (ri,„,zi,„) (si,yi) because (si,n,yi,n) (si,2/i) and /Ci'^ has a 
closed graph with IC'^'^{si) = /C(si) (see Lemma [4£j). In particular, Ti^n > 
for large n. Since is externally tangent to /C^'^ at (si,rn ?/i,n)) Lemma 
14.71 states that one can find a C^'^ tube externally tangent to /Ci n 
([0,Ti,„] X 1?^) at (n ,«) ^i,n) with 

s(n,„)=s(«,„)+(.i,„-,i,„). 

Moreover Ei^„ is also in a neighbourhood of (ri^n, -zi.n)- Lemma l48l given 
below, states that the tube /Ci H ([0, ri^„] x 1?^) is still a subsolution to the 
evolution equation for hx^ ■ Therefore, 

,71 (,'^1 .71 J 

Smce Si,„(si,„) = Si,„(ri,„) - - we have f^^i^;, ' ' ' = v^^ '^ 

and Hyl':^'^'''^ = H^l::^^^'"\ Therefore 

^Kr,yi,„)+^^"'' < F(i.J^-(^^-\i/J^;"(^^-VAi/i(yi,n,Si,„(si,„)-(zi,„-yi,„)^ . 
Letting n — )• +oo we get 

Vlly^^ + ee-^^<F{vl^^^^\X,) + \Myi,^i{si)) . (48) 
The difference between (07|) and (08|) gives 

> F(z^^2^^'^\X2) + A2/i(y2,S2(s2)) 

> {X2{1 - 9ef - Xi)h{yi,^{t) - eiu^) + ee-'\ (49) 

because, on the one hand, y.^-^ = ^^2y2)^ ^yi^^^^ = '^y2^'^^^ and Xi < X2 
(from (j44l 145145]) ) and F is elliptic and because, on the another hand, from 
[T]) and the definition of S2, 

h{y2,^2{s2)) = Hy2,^i{si) + y2-yi) > {l - 6\y2 - yi\f Hyi,^i{si)), 
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where |?/2 — yi| ^ '^^ < e. Since h > and A2(l — Oe)^ > Xi from Lemma 
14.51 there is a contradiction in ()49p and the proof is complete. 



4.3 Proofs of Lemmas \4M EM SZl and [478] 

Proof of Lemma 14.51 The proof of this lemma is close to the one of [TJ 
Proposition 4.2]. We provide it for reader's convenience. 
The assertion (i) is obvious. 

To prove (ii), we first note that, since /Ci(0) n /C2(0) = and /Ci and IC2 
are closed, there exists r > such that /Ci(t) n/C2(s) = for all < s, t < r. 
Therefore 

/i := inf{|y - x\ : < s,t < t, x G /Ci(t), y G ^(s)} > 0. 

Let e < fj, and cr < r/(2e). Then, for s G [0,r/2] and y G /Ci(s) and for any 
(t,x) G /C2, we have 



Hence dj^{s,y) > e, which proves that T"^'"" > r/2 > e. 

We prove (iii). From the definition of pa{T'^'")-, there exists t„ | 2^^'°^ 
with Pa{tn) < ee~*". Therefore there exists G K.i{tn) such that Paitn) = 
d^{tn-,yn) and, up to extract a subsequence, we can assume that yn ^ V ^ 
/Ci(T^''^) (since /Ci is closed). It follows 



Thus, we obtain the inequality pa-{T^'^) < ee"-^''"^ (note that we prove by the 
way that p^ is a lower-semicontinuous function). It remains to prove that 
the equality holds. If not, there exists y G ICilT'^''^) such that d'^CT'^'^ ,y) < 

gg-T^'"^ From the left lower-semicontinuity of the subsolution /Ci, for all 
sequence tn t {T'^'")^ ■, there exists yn G K.i{tn) which converges to y. It 
follows that 



We get a contradiction with the definition of T'^''^ and conclude for the proof 
of (ifi). 



QED 




Paitn) = rf^(tn,yn) ^ dp(r^'",y) > PaiT'^- 



Paitn) < d'^{tn,yn) < £6 



at least for n large. 
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We turn to the proof of (iv). For this we fix e > as in the proof of (ii) 
and we note that T^'" is noncreasing with respect to a. Since T*^'*^ < T*, 
limo._j.Q+ T'''" exists and we denote it by s. Note that s > e. 

We now argue by contradiction, assuming that there is a sequence (t„ — )• 

0+, yi,n e /Ci(r^''^"), (S2,n,y2,n) G ^2 SUch that 

with either S2,n = or = y2,n- Since /Ci and /C2 have closed graphs, up 
to extract subsequences, there exist yi G Ki{s) and (s2,?/2) G /C2 such that 
yi,n yi and (S2,n,2/2,n) (s2,y2)- FYom the inequahty 

e> (— (T^'"" -S2,n)' + |2/l,n-y2, 

we deduce that T^'*^" - S2,n 0. But T^''^" s > e, and so S2,n > for n 
sufficiently large, which implies that yi „ = y2,n and thus yi = 2/2- 

We now use the left lower-semicontinuity of fCi and /C2: Let t„ — )• s~. 
Since j/i = 2/2 G ^i(s) H /C2(s), there exist G A^i(in) and X2,n £ ^2 (in) 
which converge to yi = 2/2- Then 

as soon as n is sufficiently large. This is in contradiction with the definition 
of T^''^". 

QED 

Proof of Lemma 14.61 The set /C^'*^ is closed because so is /Ci. Let y G 

^i'^(si)- There is some (r, z) G /Ci such that t < si and |(r, z) — (si,y)|o- < 
e(e-^ - e-^i). Then 

^(r - < |(t, z) - (si, < e\e-^ - e-^^f < e^r - s,f . 

Since ecr < 1 (from Lemma 14.50 . we have t = si and thus y = z. So we 
have proved that /C^''^(si) C /Ci(si). The other inclusion being obvious, the 
equality holds. 

Let (s,y) G /C^''^. Then there exists (r, z) G /Ci with r < si and |(s,y) — 
< e(e~'^ — e~*M. From the definition of and Si, we have (t, z) > 

PaiT) > ee""^. It follows 

d''^{s,y) > d''^{T,z) - \{s,y) - {t,z)\^ > ee'^i. (50) 
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Taking the infimum over {s,y) G /Ci'^, we get e(/C^''^, /C2) > ee~^'^. 

Let us prove the opposite inequahty. From Lemma 14.51 we can choose 
y S /Ci(si) such that d^{si,y) = ee~'^^. But /Ci(si) C /C^''^(si). Therefore 

e(rl•^^2)<(^^(sl,y)iee-^^ 

To prove the second assertion, let yi G /Ci(si) and (52,2/2) £ be such 
that ()42p in Lemma 14.51 holds. Then obviously ()43p also holds. For such 
points, let (r, z) G /Ci be such that t < si and 

|(si,yi)-(T,z)U<e(e---e~^^). (51) 

If r < si, then, by definition of si, we have Pa{T) > ee~'^ . From the com- 
putation (l50l) with (s,y) := (si,yi), it follows d^(si,yi) > ee^"^ > ee"**! 
which is a contradiction with (|43|) . Hence r = si and therefore (|51|) gives 
(siiyi) = (ttz) G /Ci. We conclude by Lemma 14.51 

QED 

Proof of Lemma 14.71 Let be the signed distance to 9S: 

"^^^^'"^-1 -d-9Ar,z) if (M)GS. (^^^ 

Since S is of class C^'^, one can find r/ > such that is of class C^'^ in 
{|dg| < r]}, with Z?^.dg / in this set. 
Let us define S by 

E = {(r, z) G li+ X : dg((r, z) + {s, y) - (f , z)) < -e(e-^ - e"^)}. 

We first show that the tube S is externally tangent to /Ci H ([0,f] x 1?^) 
at ('r,^). At first, if (r, z) G /Ci with r < f, then by definition of JC^''^ , 
we have B(j{{t, z),e{e~'^ — e"'*^)) C /Ci"^, where 5^ is the usual open ball 
related to the norm | • l^-. In particular, since /C^''^ C S, we have d^(r, 2) < 
—e{e~'^ — e"**^). Therefore 

d-((T,z) + (s,y)-(r,z)) 

<dg(r,z) + |(s,y)-(r,z)U 

< -e(e-^ - e-"i) + e(e-^ - e''^) 

< -e(e-^ - e"^) . 

So we have proved that /Ci H ([0, f] x M^) C S. Moreover, we obviously 
have that (r,z) G dT, because (s, y) G 
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Let us show that S is a regular tube in an open interval containing r. 
For this, recall that bd($]) is defined by (j22p . If (r, z) belongs to bd(S) with 
€\e~'^ — e~'^\ < 7], then 

|d5((r, z) + (s, y) - (r, z))\ < eje"" - e"^! < r? 

and so is differentiable in a neighbourhood of (r, z) + (s, y) — (f, z) with 
-D^d^ 7^ 0. Thus S is as smooth as S. Finally, 

^(-^-)- |z).d^(.-,y-)| ^^(^-5) + ^^ ' 

since |Z)^d-(s,y)| < |Ddg(s,y)| = 1. 

QED 

Lemma 4.8 If fCi is a subsolution to the evolution equation for h\^, then, 
for any t > 0, /Ci n ([0, t] x R^) is also a subsolution for h\^ . 

Proof of Lemma \AM Let us set £i = /Ci n ([0,t] x R^). It is clear 
that /Ci is a left lower-semicontinuous tube because so is /Ci. Let S be a 
tube defined on some open time-intervall I and which is external tangent to 
/Ci at some point (s, y) with s > 0. If s < t, then, assuming without loss 
of generality that / C (0, t), S is also externally tangent to fCi and thus 

^£,)</iAi(y,s(s)). ^ 

We now suppose that s = t and, without loss of generality, that 9Sn/Ci = 
Let d^ be the signed distance to defined by ([52]) . Note that, 
since 9S n /Ci = {(s, y)}, d^ has a strict maximum on /Ci at (s, y) (at least 
on the interval /). For 7 > we introduce the mapping 

^P^{t,z) = dg(r,z) +7log(t - r) 

of class for r € /H (0,t) and when |dg(T, 2;)| is small. From standard 
arguments (see for instance the proof of Lemma 4.2 of [8j), ip.^ has a maxi- 
mum on /Ci at some point (57,2/7) G /Ci which converge to (s,y) as 7 — >• 0^ 
and such that S7 < t. Moreover, the set 

S7 = {(r, z) e (/ n (0, t)) X R^ : (^^(t, z) < ip^{s^, y^)} 

is a tube on some open interval C I (1 {0,t) with s-y G 1^, and ^7(57) 
converges in the sense to S(s). Now we note that S7 is externally tangent 
to /Ci at (57,2/7) and thus 

< /lAi (2/7' ^7(^7)) ' 
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with 
Hence 

which gives the desired inequality as 7 — )• 0"*": V^y-j < hx-i^{y,'E{s)) . 

QED 

5 Convergence to equilibria 

In this section we investigate the asymptotic behavior of the solutions of 
our front propagation problem ([6]). More precisely, we show, under suitable 
assumptions on the source S and on F, that the solution converges as i — )• 
+00 to the (weak) solution of the free boundary value problem: 

Find a set -fC G P such that hxix, K) = for all x £ dK (53) 

(recall that D is defined by (jlip ). This problem is a generalization of the 
Bernoulli exterior free boundary problem we recalled in the introduction. 
Let us first introduce a notion of weak solution: 

Definition 5.1 We say that the set K £ D is a viscosity suhsolution (re- 
spectively supersolution, solution) of the free boundary problem (in short 
FBP) 153)) for h\ if the constant tube IC{t) = K for all t > is a subsolu- 
tion (respectively supersolution, solution) of the FPP ([S]) for h\. 

Remark 5.1 There are many other definitions of weak solutions for such 
FBP: see for instance the survey paper [15|. The one we introduce here is 
the more suitable to our purpose. The idea of using sub- and supersolution 
in this framework comes back to Beurling [5]. 

In order to ensure that the FBP (I53p has a solution, we assume in the 
sequel the following: 

VA > 0, 3i? > such that Vr >R,\/x& B{0,r), hx{x,B{0,r)) < 0. (54) 

This assumption states that B{0,r) is a strict classical supersolution of the 
free boundary problem for hx for r sufficiently large. It is in particular 
fulfilled (i) when F^VjA) = Tr{A) + -Fi(z^), where Fi^v) < 0, or (ii) when 
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F = F{u) < for any v with \v\ = 1, because of the behavior of h for large 
balls (see Lemma l2.2p . Note also that the assumption implies that, for any 
ball B, 

since F is elliptic and h > 0. 

Proposition 5.2 We assume that /ilO\) and hold. Then, for any A > 0, 
there is a largest and a smallest solution of the FBP \5^) for h\, the largest 
being closed and containing any suhsolution for hx, while the smallest is 
open and is contained in all the supers olutions. 

Proof of Proposition 15.21 The proof can be achieved by a direct applica- 
tion of Perron's method. Existence and bounds for sub- and supersolutions 
are ensured by assumption (I54p and by Lemma 15.31 below. 

QED 

Lemma 5.3 We assume that Iil0\) and |5^[ j hold. Let X > be fixed. Then 
there exist e > and R > such that, if K is a subsolution ( respectively 
a supersolution) of the FBP ^5^^ for h\, then K C B{0,R) (respectively 
Se C K where is defined by (14^)- 

Proof of Lemma 15.31 Let r be the smallest nonnegative real such that 
K C B{0,r). Then there is some point x E dK D d{B{0,r)). Using the 
constant tube -6(0, r) x 12 as a test-tube and the fact that K is a subsolution, 
we have hx{x, B{0,r)) > 0, which implies that r < R where R is given by 
Therefore we have K C B{0,R). 
The assertion for the supersolution can be proved in a similar way, by 
using assumption (llOp and Lemma Essaying that h{x, S'e) is large for e > 
small and x G dS^. 

QED 

Next we address the uniqueness problem. The main assumption for this 
is that S is starshaped. We also suppose that F = F{v, A) satisfies the 
subhomogeneity condition: 

F{v,-iA)>-iF{v,A) V7>1, (55) 

and that the following compatibility condition between F and S holds: 

F(i/f,Ff)<0 Vxg55. (56) 
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Assumption ()55p is fulfilled for instance if (i) F{i',A) = Tr{A) + Fi{i'), 
where Fi{i') < 0, or if (ii) F = F[v) < for any with = 1, while 
assumption (j56]) is always satisfied for F as in (ii), and is satisfied for sets 
with negative mean curvature if F is as in (i). 

Theorem 5.4 Let us assume that S is strictly starshaped at 0, with G 
int(S'), g = 1 and that ilO\} . ^54\ ), i55\) and i56\) hold. Then, for any A > 0, 
the FBP ( 153)) for h\ has a unique solution denoted K\. Moreover, K\ is 
starshaped at for any A > and the map A — )• Kx is continuous for the 
Hausdorff topology. 

Remark 5.2 

1. Uniqueness of solution means that, if Ki and K2 are two solutions of the 
FBP for hx, then = and ini{Ki) = int(K2). 

2. Such a uniqueness result is classical in the literature, see in particular 
Beurling [5], Tepper [27] and the survey paper |15|. 

In order to prove Theorem [23] we need three Lemmas. The first one explains 
that the homothetic of a subsolution is still a subsolution. The second one 
allows to compare sub and supersolutions of the FBP. The last one shows 
that subsolutions of the FBP for hx when A is small are necessarily close 
to S. 

Lemma 5.5 Assume that S is starshaped at 0, g = 1 and that 155\) holds. 
If K is a subsolution of the FBP Ii53\) for hx, then pK is a subsolution of 
the FBP for hpx for any p € (0, 1) such that S CC pK. 

Proof of Lemma 15.51 For sake of clarity, we do the proof in a formal 
way, by assuming that K is smooth. If not, it is enough to do the same 
computations for the test-surfaces. We first notice that 

h(x, pK) > \h{-,K) yp G (0, 1), Vx G d{pK) . (57) 
P 

Indeed, if u is the solution to ([9]) with K instead of then v{x) = u{x/p) 
is a subsolution of equation ([9]) with pK instead of Q (we use here the fact 
that S is starshaped, that g = 1 and, thus, that < u < 1). Then 

h{x, pK) > \Dv{x)\'^ = ^\Du{x/p)\^ = ^h{x/p, K) Vx G d{pK) . 
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Next we also notice that = v^j^ and Hx^ = jH^j^. Hence, for any 
X G d{pK), we have 

hpx{x,pK) = F{uf,HP^) + p\h{x,pK) 

> ni^^/,^j;H^/,) + fhix/p,K) 

> il/p)hxix/p,K)>0 

because K is a subsolution for hx- Hence pK is a subsolution for hpx. 

QED 

Lemma 5.6 We assume that jj^p holds. Let 0<A<A, i?>0 and 7 > 
be fixed. Then, there is a constant n > 0, such that for any A < Ai < A2 < A, 
for any subsolution Ki of the FBP \5^) for hx^ and any supersolution K2 
for hx2 with 

CC K1CCK2 CC B{0,R--/), 

we have 

Ki + k{X2- Xi)B{0,l) C K2 , 

where the sum in the above inclusion denotes the Kuratowski sum between 
sets. 

Proof of Lemma 15.61 Let ^ > be the constant given by Lemma 12. 4[ 
From the assumption Ki CC K2, we have Ki n R^\K2 = and we can 
find yi £ Ki and y2 G M^\K2 such that 

\yi-y2\= mill \zi - Z2\ . 

zi<^Ki,Z2£R^\K2 

Without loss of generality we can assume that 1/2! < 1/^, since otherwise 
the result is obvious. 

Using now the interposition and approximation results (see Proposition 
13.11 and Theorem 13. 2p , the fact that Ki is a subsolution for hx^ and K2 a 
supersolution for hx2 , and proceeding as in the proof of Theorem 14.11 one 
can find an open set J: C with C^'^ boundary and (iV - 1) x (iV - 1) 
matrices Xi < X2 such that 

0< F(z.^^,Xi) + Ai/i(yi,S) (58) 

and 

^>F{vl+y^-y\X2) + X2h{y2,^ + y2-yi) • (59) 
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Since 
using Lemma [2 



S+J/2-J/1 

^y2 



and Xi < X2, we get by subtracting ([58|) to ([59|) and 



> 



(60) 



In order to complete the proof, we have now to check that h{yi, S) is posi- 
tive. By Hopf maximum principle, we just have to show that the connected 
component T,' of E which contains yi has a non empty intersection with the 
source S. For this, we argue by contradiction, by assuming that S' n S" = 
(see Figure m for an illustration). Let K[ := Ki n S'. Note that 




e{K[,Y\T.') > e(S', S\S') > (61) 

since S is bounded with a smooth boundary. Let z £ S and y[ be a point 
of maximum of the euclidean norm | • —z\ on K[. Note that \y[ — z\ > 0. 
The ball B := B{z,\y[ — z\) is externally tangent to K[ at y[. Thanks to 
(f6T]) . one can build an open set T,[ with a boundary, such that K[ C 
T,[ n = ^, y[ £ dTi'^, and for which there is a neighboourhood O of y[ 
with T,[r\0 = B no. Note that S[ n S" = since 5 CC E\E'. Let us set 
Si = S'^U (S\S'). Note that Si is externally tangent to Ki at y[. Moreover 

/iA(y'i,Si) = F(z.J/,i7j/) + A%l,Si) = F{u^^,H^^) 

since /i(yi,Si) = 0. By (f54j) . for all r > i? and all x G dB{0,r), we have 
/iA(a;,5(0,r)) < 0. Therefore F(i/f //f ^°''^^) < since h{x,B{0,r)) > 0. 
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By ellipticity, we have F{i>^,^^'^ \h^,^^'^ ^) < even for r' < R and \x'\ = r' 
since u^,^^''' ^ = f^^M- ^* follows 

/iA(yl,Si) = F(z^g,i?,^)<0 

which is a contradiction since Ki is a subsolution. So > 0. 

Then ([60]) leads to inequality 

\yi - y2\ > k(A2 - Ai) 
where k := 1/(20 A), which completes the proof. 

QED 

Lemma 5.7 Under the assumptions of Theorem \5.4\ for any e > 0, there 
is Xq > such that, for any subsolution K of the FBP ( 153|) for hx with 
A e (0, Ao), we have K C (see (T4\ ) for a definition of S^). 

Proof of Lemma 15.71 From assumption ()56p and the regularity of the 
boundary of S, there is some a > such that 

F{iy^,H^)<-a e dS . 

Let us notice that a similar inequality also holds for pS, for p > 1, because 

thanks to assumption (fSSl) . 

Let us now fix e > and po > 1 such that p^S C S^- Note that 
S CC pqS since S is strictly starshaped. We set Aq = a/(/3po)i where 
j3 = sup^.g9(pp5) h{x,poS). Let if be a subsolution for hx with A G (0, Aq). 
We denote by p > 1 the smallest real such that K C pS. In order to prove 
that p < po, we argue by contradiction and assume that p > po. Since pS is 
externally tangent to K at some point x G dK and K is a subsolution, we 
have 

0</..(...5)^F(.f,«f) + AM...5)<-2 + AM..pS) 
where, from inequality ([57|) . 

Hence < — ^ + A (^^) which implies that p < AoPo/^/a = po, a contra- 
diction. So we have proved that p < pQ. Therefore K C poS C Sf.. 
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QED 



Proof of Theorem 15.41 Let us denote for any A > by Kx the maximal 
solution of the BFP for hx- Note that A — )• Kx is nondecreasing since Kx 
contains all the subsolutions for hx- 

We first check that Kx is starshaped at 0. Indeed, from Lemma |5.5|. for 
any p £ (0, 1) sufficiently close to 1, the set pKx is a subsolution for hpx, 
and thus for hx- Since Kx contains all the subsolutions, we have pKx C Kx 
for any p £ (0, 1) sufficiently close to 1, which implies that Kx is starshaped. 

Next, we show that the map A — t- Kx is continuous for the Hausdorff 
topology. From the stability property of solutions (Proposition 14. 4p . the 
decreasing limit of the Ky converges to Kx when A' — )• A+. Hence we 
only have to show that Limx/^x- ^X' equals Kx, where Lim denotes the 
Kuratowski limit (see (|39|) ). 

Since, for any p £ (0, 1) sufficiently close to 1, the set pKx is a subsolution 
for hpx, we have pKx C Kpx- Therefore 



So we have checked that A — )• Kx is continuous. 

Let us finally prove that, for any A > 0, Kx is the unique solution of for 
hx- Let K be another solution. Note that K C Kx- From Lemma 15.71 we 
can find some Ai > such that Kj CC K because S CC K- Let us set 



We now use Lemma 15.61 with r > and R such that Sr C Kj_^ and Kx C 
-6(0, R)- There is a constant k > such that for any Ai < A' < A, 

Ky + k(A - X')B C K . 

The continuity of the map A' — t- Kx' then implies that 

Kj + k{X-X)B cK. 

Therefore A = A since, otherwise, the continuity of A' — )• Ky would also 
imply the existence of e > such that Kj^^ CC K, a contradiction with 
the definition of A. Therefore Kx C K- 

In order to prove that int(i^) = mt{Kx), we notice that mt{Kx) = 
Ua'<a ^A') because Kx' CC Kx for A' < A, and therefore the equality A = A 
implies that mt{Kx) C int(i^'). Since the converse inequality is obvious, the 
proof of Theorem 15.41 is complete. 



Kx = Lim pKx C Lim Ky C Kx- 



A = sup{A' \Ky CCK} . 
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QED 



Corollary 5.8 (Asymptotic behavior) Under the assumptions of Theo- 
rem \5.4\ if fC is a solution of the FPP ^ for h\, then IC(t) converges, for 
the Hausdorff metric as t ^ +00, to the unique solution K\ of the FBP 
\5^) for h\ while K{t) converges to R^\K . 

Remark 5.3 

1. Note that the above result holds for any solution lC{t) of the FPP ^ 
with any initial position /C(0) G T>. 

2. The proof of the asymptotic behavior which follows relies strongly on the 
uniqueness of the solution of the limit problem (f53l) . 

Proof of Corollary 15.81 Let us fix Ai < A < A2. From lemma [2T3l and 
there are r > and i? > such that 5,. and B{0,R) are respectively 
sub- and supersolution to the FBP for hx-^ and hx.2- We can also choose 
r > sufficiently small and R > sufficiently large so that 5,. CC /C(0) CC 
B{0,R). The inclusion principle then states that 

Sr CCJC{t) CC B{0,R) Vt>0. 

Let K* be the Kuratowski upperlimit of /C(t) as t — t- +00 (see ()39p ). Note 
that Sr C K* C B(0,R) and that the constant tube R x K* is actually the 
upperlimit of the solutions /C(- + r) as r — )• +00. From the stability of 
solutions (see Proposition I4.4|) . the constant tube M x K* is a subsolution 
of the FPP for hx. Hence, K* is a subsolution of the FBP (|53p for hx and 
we have K* C Kx- 

In the same way, if we set L* to be the Kuratowski upperlimit of K{t) 
as t +00, then R^\L* is a supersolution to FBP ([53]) for hx- Since Kx 
is the unique solution for hx, Kx is also the smallest solution, which implies 
that Kx C R^\L*. Hence 

K* C KxC W^L*. 
Since R^\L* C K*, the proof is complete. 

QED 
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